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The charmed-strange meson mases are calculated on a quenched lattice QCD. The charm and 
strange quark propagators are calculated on the same lattice with the overlap fermion. 16^ x 72 
lattices with Wilson gauge action at j3 = 6.3345 are used. The charm and strange quark masses are 
determined by fitting the J/y/ and (j) masses respectively. The charmed strange meson spectrum 
for the scalar, axial, pseudoscalar and vector channels are calculated. They agree with experi- 
ments. In particular, we find the scalar meson mass to be 2248(78) MeV which is in agreement 
with that of £);o(2317). 
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In 2003, BaBar Collaboration announced the discovery of a charmed strange meson £)*q(2317) 
[|l|]. CLEO also reported the observation of this particle in the same year In the following 
plot, we show the masses of the charmed strange mesons from the Particle Data Group, in which 
the newly discovered D*q(2317) is a scalar meson and Dsi(2460) and Dii(2536) are axial-vector 
mesons. 



Figure 1: The charmed strange meson spectrum from PDG 
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Predictions of these charmed strange meson spectrum have been made in the quark model |^]. 
While they gave good prediction for the tensor and axial-vector, pseudoscalar, and vector 
mesons, its prediction of the 'Pi axial-vector at 2.53 GeV is '-^ 70 MeV above the experimen- 
tal Dvi(2460). More puzzling is the prediction of the scalar meson at 2.48 GeV which is ~ 160 
MeV above D*q(2317). This discrepancy has prompted the speculations that D*q(2317) is a DK 
molecule [Q] or four quark state or a threshold effect instead of a cs meson. 

There are also a few lattice calculations. Lattice NRQCD calculation with quenched approxi- 
mation gives m{D*Q) = 2.44(5)GeV [^]. The Uf = 2 calculation with the heavy quark at the static 
limit gives m{D*Q) = 2.57(ll)GeV [^. These are also significantly heavier than the experimental 
mass of 2.317GeV. The recent calculation with a relativistic heavy quark ( RHQ ) action gives 



Am = m(D*)) -m(D,) = 0.1243(28)GeV, or m{D*Q) = 2.093(3)GeV [|TO|l, which is significantly 
lower than the experimental mass of D*q. 

Although, in principle, lattice QCD is an ideal tool to calculate hadron spectrum from the first 
principle, in practice it suffers from systematic errors such as due to discretization effects. These 
errors can be large for fermion actions which do not have chiral symmetry at finite lattice spacing. 
In particular, the ma errors can be substantial for heavy quarks with the commonly used lattice 
spacing a ~ 0.1 fm. 

The overlap fermion action obeys chiral symmetry at finite lattice spacing and is, thus, free 
of 0{a) and 0{ma) errors. It is shown that the effective quark propagator of the massive overlap 
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fermion has the same form as that of the continuum [11], i.e. 

/I 

Dov{m) = Dov+ma{\ - -Do 
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Dov{m) 
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where Dc 



1-4 A 



and{D„75} = 0. 



As it can be seen from Eq. the formalism is the same for all quark masses. As such, it 
can be used for both light and heavy quarks as long as the 0{m^a^) errors are negligible [11|. In 



fact, the 0{a^) error is also small [12|. By examining the dispersion relations and the hyperfine 
splittings, it is found that one could use ma < 0.5 and still keeps the 0{m^a^) errors to less than 



3% to 4% [13|. 



We carried out a study of charmonium on a quenched 16^^ x 72 lattice with Wilson gauge action 
at j8 = 6.3345. With the ro = 0.5fm scale we obtain a = 0.0560fm. Multi-mass overlap inverter is 



used to calculate propagators for 26 quark masses for ma from 0.02 to 0.85 [ 14] 



Figure 2: Quark masses used to calculate // 1// and ^ masses on the 16x12 lattice. 
16^ X 72 with Overlap Fermions 
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Matching the J / Xjf mass m{J/ Xjf) = 3.097GeV, we obtain the charm quark mass in lattice unit 
to be nica = 0.431(1). Inputting the (^>(1~) particle mass m{^) = 1.020GeV, we obtain the strange 
quark mass in lattice unit to be niga = 0.0205(32). The vector meson mass for a range of quark 
masses is shown in Fig. 2. We find that the charm and strange quark masses are in the range 
0.01 < MqU < 0.5. Thus the 0{m^a^) error is expected to be around 1% for the masses of the 
charmed strange mesons.. 

We construct the charmed strange meson correlators with the standard local interpreting fields: 



0" ^ X{x) = i\f{x)ysW{x) (3) 

0+ ^ X{x) = W{x)W{x) (4) 

I- ^X{x) = xi/{x)Yj Wix) y = 1 , 2, 3 (5) 

1+ =^ Xa{x) = xi/ix)Y5YMx) j = 1,2,3 (6) 

1+ ^ Xb{x) = n^nrMx) {ij} = {12}, {23}, {31} (7) 



We do not consider the mixing between Xa and Xb for the axial- vector mesons in this work. 

All the meson correlators with 100 configurations show reasonably nice cosh behavior, such 
as in the scalar channel in Fig. 3. 

Figure 3: The charmed strange meson correlator in the scalar channel. 
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We used 5 quark masses each around those of the charm and strange to construct 25 correlators 
for the above interpolation fields. We first fix the light quark mass and vary the heavy quark masses 
to fit the meson masses which correspond to the charm mass at nica = 0.431. The results for the 
scalar and vector mesons are plotted against the Ught quark mass in Fig.4. 

We then fix the heavy quark at the charm mass and fit the light quark masses which correspond 
to the strange mass at m^a = 0.0205. This way, we obtained the charmed strange meson masses for 
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Figure 4: The charmed scalar and vector meson mass varying with light quark mass 
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each of the meson interpolation field. The numerical results are shown in Table 1. The experimental 
data are from the PDG particle hstings. 

Our lattice results are plotted in Fig. 5 together with the experimental data. 
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Table 1: charmed strange meson masses 



Particle 


Mass( xa ) 


Lattice (MeV) 


Exp. (MeV) 


Ds{0-) 


0.5608(31) 


1976(11) 


1968.49(34) 


D*in 


0.6049(36) 


2131(13) 


2112.3(5) 


£».:o(o+) 


0.638(22) 


2248(78) 


2317.8(6) 


DM) 


0.684(18) 


2410(63) 


2459.6(6) 




0.703(26) 


2476(92) 


2535.35(84) 



Figure 5: Calculated charmed strange meson spectrum in comparison with the experimental results 
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We see that the pseudoscalar (Dj) and the vector (D*) agree with the experimental results very 
well, assuming the D* at 2112 MeV is the 1~ meson. All the p-waves states (0"*" and the two 
1+) agree with the experiments within one sigma. In particular, we find the 0+ meson with the 
Ij/Xj/ interpolation field at 2248(78) MeV to be much higher than the quark model prediction. This 
suggests that D*q(2317) may well be the conventional cj meson. 

To have a better estimate of the D*q mass, we used the ratio method to calculate the mass dif- 
ference. Since the correlators are calculated from the same gauge configurations, the ratio method 
is expected to reduce fluctuations. We calculated the ratio between the scalar and pseudoscalar 
correlators. 



G(0+) 
G(O-) 



exp(—Am-t), Am = m(0^)— m(0 



(8) 



The preliminary result is Am = m{D'^) —m{Ds) = 275(53) MeV which is about a sigma from the 
experimental result of 349.3(7) MeV. 

To summarize, we have carried out a calculation of the charmed strange meson masses with 
the overlap fermion on quenched 16^ x 72 lattices. The calculated masses are consistent with the 
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experimental results within error bars. The = 1~ meson mass matches with D*(?)(2112) very 
well. This impMes that D*(?)(2112) is indeed the expected vector meson. The scalar meson mass 
at 2248(78) MeV is higher than the prediction of the quark model and is in agreement with the 
experimental mass of £)*q(2317). 

This work is partially supported by DOE Grands DE-FG05-84ER40154. We thank H.Y. Cheng 
for bringing the problem of Z)*q(2317) to our attention. 
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